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ABSTRACT 

i— h \ In equivariantly gauge-fixed SU (N) Yang-Mills theories, the gauge 

p symmetry is only partially fixed, leaving a subgroup H C SU (N) un- 

fixed. Such theories avoid Neuberger's nogo theorem if the subgroup H 
contains at least the Cartan subgroup [/(l)^ -1 , and they are thus non- 
perturbatively well defined if regulated on a finite lattice. We calculate 
| the one-loop beta function for the coupling g 2 = £g 2 , where g is the gauge 

coupling and £ is the gauge parameter, for a class of subgroups including 
the cases that H = [/(l)^ 1 or H = SU{M) x SU(N - M) x 17(1). The 
coupling g represents the strength of the interaction of the gauge degrees 
of freedom associated with the coset SU(N)/H. We find that g, like g, 
is asymptotically free. We solve the renormalization-group equations for 
the running of the couplings g and g, and find that dimensional transmu- 
tation takes place also for the coupling g, generating a scale A which can 
be larger than or equal to the scale A associated with the gauge coupling 
g, but not smaller. We speculate on the possible implications of these 
results. 



*Permanent address: Department of Physics and Astronomy, San Francisco State University, San 
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1. Introduction 



In a recent paper pQ we constructed equivariantly gauge-fixed SU(N) Yang-Mills 
(YM) theories on the lattice, generalizing earlier work for the gauge group SU(2) [2|. 
The aim was to gauge fix as much of the SU(N) gauge symmetry as possible while 
avoiding a nogo theorem stating that the standard BRST procedure breaks down 
non-perturbatively on the lattice in the compact formulation jHj. It turns out that it 
is possible to gauge fix SU (N) down to the maximal abelian subgroup H = ^/(l)^ 1 
while avoiding the theorem. 

There are other ways to gauge fix YM theories non-perturbatively, but usually the 
gauge-fixing procedure is extraneous to the path integral (for a review, see ref. 
Our aim was to maintain explicit locality as well as a version of BRST invariance. 
At the price of leaving (at least) the Cartan subgroup unfixed, equivariantly gauge- 
fixed theories fulfill this aim: Their complete definition consists of a path integral 
with a local action which, in turn, is invariant under an equivariant BRST (eBRST) 
symmetry. They are rigorously well-defined on a finite lattice, and one can prove that 
correlation functions of gauge-invariant operators in the un-fixed theory are exactly 
the same as those in the equivariantly gauge-fixed theory. The eBRST symmetry 
is still nilpotent in the following sense: The eBRST algebra closes on the subgroup 
H, which is left unfixed. This implies that on if-invariant operators eBRST is still 
nilpotent, and thus the usual consequences of nilpotency still apply. 

Equivariantly gauge-fixed SU(N) YM theories, in which a subgroup H containing 
at least the Cartan subgroup of SU (N) remains unfixed, are thus non-perturbatively 
well defined, local gauge theories, possessing an exact BRST-like invariance. Of 
course, if one would like to completely gauge fix the theory non-perturbatively, the 
remaining subgroup H would still need to be fixed. At present, it is not known how 
to do this such that both the nogo theorem is avoided and a form of BRST symmetry 
is preserved. 1 In contrast, in perturbation theory, complete gauge fixing is always 
possible, and, in fact, necessary. In the equivariant case, we only need to further fix 
the subgroup H in order to set up perturbation theory. This leads to modifications 
of the equivariant Slavnov-Taylor identities, but these can be taken into account. 
Unitarity and renormalizability to all orders are still guaranteed (for the case that H 
is the Cartan subgroup, this was investigated in detail in ref. [T]). 

In this paper, we want to start a dynamical investigation of the phase diagram of 
equivariantly gauge-fixed YM theories (without any fermions or other matter fields). 
We would like to know whether any phases exist other than the one containing the 
continuum limit corresponding to the standard YM universality class. In particular, 
we would like to know whether any phases exist in which eBRST invariance and/or 
any of the other global symmetries exhibit spontaneous symmetry breaking (SSB). 
We emphasize that this is a completely non-perturbative issue: equivariantly gauge- 
fixed YM theories reproduce standard perturbation theory JX (including perturbative 
unitarity), and the perturbative expressions for any gauge- invariant quantity are thus 
equal to those obtained in YM theory in any standard perturbative gauge. 

^or a way out for purely abelian theories, see ref. It is not known how to generalize this 
trick to our case, where the abelian degrees of freedom are only part of the compact gauge field 
represented by the link variables. 
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Let us first assess the fact that also non-perturbatively the standard YM uni- 
versality class can be obtained from its equivariantly gauge-fixed cousin. As already 
stated above, gauge-invariant correlation functions remain unchanged after equivari- 
ant gauge fixing on the lattice in a finite volume. In this situation, it is therefore clear 
that whatever happens in the gauge- fixing ( "unphysical" ) sector of the theory cannot 
"spill over" to the physical sector: the physics should be that of the unfixed theory. 
If this is true in any finite volume, it remains true in the infinite- volume limit, and 
we therefore recover the standard YM universality class. 

This raises the question whether SSB could occur at all in such a way as to modify 
not only the gauge- variant, but also gauge-invariant quantities. As is well-known, in 
order to investigate possible SSB in any finite volume, an external field needs to be 
turned on in order to "seed" the SSB, and this external field is only taken away after 
the infinite- volume limit is taken. Since this external field will break eBSRT symmetry 
explicitly, it is not a priori clear what the dynamical effects of SSB will be, even for 
the gauge invariant sector of the theory. One anticipates that the parameters of the 
theory (the gauge coupling g and the gauge- fixing parameter £) can be chosen such 
that the usual YM universality class can be recovered. However, whether this will 
happen or not is a question for the dynamics of the theory: there is no obvious reason 
that the vanishing-external-field and infinite-volume limits will commute everywhere 
in the phase diagram. A separate question is whether any phase with SSB might or 
might not contain any critical point corresponding to an interesting continuum limit. 

We take as our primary point of view that it is interesting to understand the 
phase diagram of an equivariantly gauge-fixed YM theory in its own right. This 
is corroborated by our main result, that also the coupling constant governing the 
longitudinal sector of our theory is asymptotically free. Asymptotically-free couplings 
are rare, and suggest strong-coupling dynamics, like bound states and/or SSB, at low 
energies. 

Another motivation for our interest in the phase diagram of equivariantly gauge- 
fixed YM theories comes from our interest in the non-perturbative study of chiral 
gauge theories with the help of the lattice [S] . Lattice- regulated chiral fermions break 
the gauge (and, hence, any BRST-type) symmetry explicitly [II El II], and it is there- 
fore important to explore the possible dynamics of equivariantly gauge-fixed theories 
in a wider context than that in which the infinite-volume limit is taken while main- 
taining eBRST symmetry rigorously in the regularized theory. 

Let us consider this in a little more detail. On the lattice, with chiral fermions, 
eBRST is broken by the regulator, i.e. by irrelevant terms (defined as terms in the 
lagrangian which vanish in the classical continuum limit). Because of this, counter 
terms are required |B1 E], whose coefficients will need to be tuned in order to 
recover a quantum continuum limit. A number of these counter terms break eBRST, 
and their coefficients can thus be considered external fields in the sense discussed 
above. The same reasoning thus applies: It is a dynamical question whether SSB of 
eBRST (and/or any other global symmetries) takes place, and what the consequences 
of this might be for gauge-invariant quantities. 

Returning to the pure YM case, the number of gauge couplings in the gauge- 
fixing sector depends on the subgroup. If we choose H = [/(l)^ -1 or H = SU(M) x 
SU(N — M) x £7(1) (which is a maximal subgroup of SU(N)), the theory has just 
two couplings, which we will take to be the square of the gauge coupling, g 2 , and the 
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product of the gauge- fixing parameter £ and g 2 , denoted by g 2 = £g 2 . We will refer to 
this latter coupling as the "longitudinal" coupling, since it determines the dynamics 
of the "reduced" theory [T] in which the gauge field is restricted to be pure gauge. 
As we already said, in order to study the phase diagram a non-perturbative study 
is required. Such a study is outside the scope of this paper. But, it is imperative, 
and interesting, to begin with an investigation of the scaling behavior of the theory 
in the weak-coupling regime (small g 2 and g 2 ). In this paper, we calculate the beta 
function for g 2 to one loop (the one-loop beta function for g 2 is known of course), 
and solve the resulting renormalizat ion-group (RG) equations. It turns out that also 
the coupling g 2 is asymptotically free, and, just like g 2 , generates its own associated 
infrared scale. 

For the calculation of one-loop beta functions, we may use the continuum theory 
because they are independent of the regulator chosen, and thus most easily calcu- 
lated using dimensional regularization. Therefore, in Sec. 2, we summarize only the 
continuum form of the theory (for the lattice version, we refer to ref. [Tj). In Sec. 3, 
we present the one-loop beta functions, and we solve the resulting RG equations in 
Sec. 4. In Sec. 5 we speculate on the interpretation of our results. We conclude in 
Sec. 6. In Appendix A we summarize some group-theoretical details, and discuss a 
wider class of subgroups. 



2. Lagrangian for equivariantly gauge-fixed SU(N) theory 

We start with the lagrangian for an equivariantly gauge-fixed SU(N) YM theory. 
The YM lagrangian for a gauge field = V^T a with gauge coupling g can be written 
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£ym = ^tr(F 2 u ), (2.1) 
iF^ = [Dp{V),D v {V)\ , 

We take the hermitian generators T a of the group to satisfy tr (T a T b ) = ^5 a t> and 
define the structure constants f a b c through [T a , T b ] = if a b c T c - The structure constants 
are completely anti-symmetric in all three indices. 

In order to define the gauge-fixing and ghost parts of the lagrangian, we divide 
the generators into two sets. One set, {T*}, generates the subgroup H, which is left 
unfixed, while the remaining generators, {T a }, span the coset G/H, with G = SU(N). 
We will use indices i, j, k, ... to indicate H generators, a, j3, 7, . . . for Gj H generators, 
and a,b,c, . . . for generators of the full group, G. We have that 

fiaj fodj fija , (2.2) 

because the product of two elements of H is again in H. 
Correspondingly, we split the gauge field into two parts, 

V, = A^ + W,, (2.3) 
A, = 4T\ W, = W«T% 



^Throughout this paper we work in euclidean space. 
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and we introduce coset-valued ghost and anti-ghost fields 



C = C a T a , C = C a T a . (2.4) 

The gauge-fixing condition FiV) has to be covariant under local H transformations. 

The simplest generalization of the usual Lorenz gauge is obtained by taking the 
derivative to be i7-covariant, and we thus choose 

F{V) = V^{A)W, = d,W, + i[A„ Wg . (2.5) 

The gauge-fixing lagrangian is then given by [I] 

A?/ = J 2 tr {V,{A)W,f (2.6) 



-2 tr (CV^A)V^A)C) + 2 tr ([W l ^~0\ H [W fl , C) H 
+i tr ((V^A)C)[W^,C} + [W fl ,C](V^(A)C) 



^ ( tr {(° 2 ) S/ n (° 2 ) g/ n) + I tr C ^nf + tr (({<7, C} H ) 

Here the subscripts TC and Qjli indicate projections onto the linear space spanned by 
the generators of H and G/H, respectively. The full lagrangian, obtained by adding 
Eq. (|2.1j) and Eq. (J2.6)) . is invariant under (on-shell) equivariant BRST (eBRST) 
transformations, local H transformations, and various other global symmetries. A 
complete discussion of these symmetries, the off-shell form, and the construction of 
this lagrangian can be found in ref. [I]. We note that if is chosen to be pure gauge, 
so that F^ v vanishes, indeed the "longitudinal" coupling g 2 governs the dynamics of 
the theory, since it is the only coupling which occurs in Eq. ()2.6j) . 

In order to derive the Feynman rules for this theory, several steps are still neces- 
sary. First, in order to develop a perturbative expansion in g, we rescale the gauge 
fields, A^ — > gA^, — > gW^. Doing this, the coupling in front of the quadratic 
term for in Eq. ()2.6|) becomes 

£ g 2 

in which we identify £ as the usual gauge-fixing parameter. The W propagator in 
momentum space is given by 

Note that the coupling g 2 still multiplies the four-ghost terms after the rescaling of 
the gauge fields. They will play no role in this paper. 3 Then, in order to define the 
A^ propagator, also the remaining gauge symmetry H needs to be fixed. We do this 
by choosing simply a straightforward Lorenz gauge for the field A^, i.e. we add 

1 9 

C H = - tr (d^Y + F-ghost terms (2.9) 



3 In perturbation theory, these couplings play a role in keeping the ghosts massless [5]. 
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(after the redefinition —>■ gA^) , where a is another gauge- fixing parameter. If H is 
non-abelian, an if -ghost sector needs to be added; in Eq. (|2.9|) the if -ghost sector is 
the standard one. However, for our one-loop calculation we will not need this sector. 
We emphasize that if gauge fixing is only done in order to set up perturbation theory; 
the non-perturbative theory we are interested in is the (lattice version of) the sum 
of Eq. (j2.1|) and Eq. (|2.6J) . without Eq. (|2.9J) . We also recall that the construction of 
equivariantly gauge-fixed SU(N) YM theories goes through on the lattice, with all 
symmetries preserved except space-time symmetries (which are reduced to the usual 
lattice space-time symmetries). For details, we refer to ref. [T]. 

3. One-loop beta functions 

In order to obtain the one-loop beta function for g 2 , we need to calculate the scale 
dependence of the gauge parameter £. This can then be combined with the known 
beta function for the gauge coupling g 2 , 

in order to obtain the beta function for g 2 . Here t = log(/i/A), with \x the renormal- 
ization scale, and with A the integration constant. 

The scale dependence of £ can be found from the ratio of the transverse and lon- 
gitudinal parts of the vacuum polarization for the coset gauge field W. The following 
types of diagrams contribute to this vacuum polarization. First, there are diagrams 
with the usual vertices from Eq. (j2.1|) . which lead to contributions proportional to 

facdf/3cd = C2^af3 = N5 a p , (3.2) 

cd 

where C2 is the quadratic Casimir for the adjoint representation, equal to N for 
SU(N). Then, there are vertices coupling the to the if gauge field A^ coming 
from the first term in Eq. (|2.6|) . Therefore, there are extra contributions involving 
only the structure constants / aj7 (recall that /y 7 = 0). The combination that appears 
in the W vacuum polarization is 

y faijfpi-y = H a 5 a /3 , (3.3) 

which defines the constants H a . Here, we will restrict ourselves to cases where Ch = 
H a is independent of a. 4 For the maximal abelian subgroup H = [/(l)^ 1 one finds 
that this is the case, with Ch = 1- Another example where H a turns out to be 
independent of a is the (maximal) subgroup H = SU(M) x SU(N — M) x £7(1), 
for which Ch = N/2. An additional contribution involving Ch comes from the third 
term in Eq. ([2.6)1 . because of the TC projections on the commutators; however, at 

4 In these cases, one may verify that there is indeed only one longitudinal coupling g 2 , meaning 
that each component of W M has the same transverse and longitudinal renormalization constants. 
For H = [/(l)^ -1 this follows from invariance under the "skewed" permutation group discussed in 
ref. pQ; for H = SU(M) x SU(N — M) x U(l) it follows from invariance under H. For some other 
subgroups, we refer to Appendix A. 
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one loop this is a tadpole that does not contribute to wave-function renormalization 
constants. Finally, there are ghost-loop contributions coming from the fourth and 
fifth terms (third line) in Eq. (J2.6|) . involving only the structure constants f ai s, in the 
combination 

^ fa-ysffi-yS = G a 8 a /3 , (3.4) 

which defines G a . The constants H a and G a are related, because 

C*2 — X/ facdfpcd = 2 ^ fod-yfpi-y + X] faySfp-yS = {^H a + G a )8 a p (3-5) 

17 7<5 

(we used again that / aii = 0, c.f. Eq. It follows that H a < N/2. If C H = H a 

is independent of a, so is Cg = G a , and C2 = 2Ch + C G with Ch < N/2. For 
# = U(l) N ^ we thus have that C G = N - 2; while C G = for H = SU{M) x 
SU(N — M) x £7(1) (all f a p 7 = 0). For some group-theoretical details, see Appendix 
A, where we prove Eq. ()3.3|) (and thus Eq. (|3.4|) ) and calculate H a for some examples. 

Since the A field does appear on the loops contributing to the W vacuum po- 
larization, we need the A propagator. This is why we fixed also the H gauge, as 
already discussed in the previous section. If H is non-abelian, H ghosts appear, but 
they do not couple to the W fields, and therefore do not contribute to the W vac- 
uum polarization at one loop. Also, we expect the one-loop beta function for g 2 to 
be independent of the second gauge parameter a, as is indeed borne out by explicit 
calculation. The one-loop beta function for g 2 is thus independent of the gauge used 
to fix the remaining gauge group H. 

The actual calculation of the W vacuum polarization is straightforward. Using di- 
mensional regularization and minimal subtraction, we find for the one-loop transverse 
and longitudinal wave-function renormalization constants 

Ztrans = 1 + 7^ ( ( y " £) C G + - £ - 2 «) Ctf) , (3.6) 




I6n 2 e \ \t Ca+ [l 



Zion 9 = l + 77^r \ — £C G + \- + 3-2a + £\C H 



The renormalization constant relating bare and renormalized couplings by £ = 
Z^ £ is equal to the ratio of Z trans and Zi ong , and combining this with Eq. (j3.1j) we 
obtain the g 2 beta function: 5 



17,2 
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bg 4 + cg 2 g 2 + df) , (3.7) 



dt 16vr 2 

b = 6C H , 

c = b-j C G -~C H = 12C H + 3C G = 3N + 6C H , 

d = 2C H + ^C G = ^N + C H . 

From this result we conclude that the coupling g 2 is asymptotically free, just like 
g 2 , because all three coefficients b, c and d are positive. Unlike Eq. (|3.1|) though, 

5 For N — 2 our result agrees with the result obtained in ref. F° r more general subgroups 

than the ones considered here, Ch and Cq in Eq. (|3.7|) get replaced by H a and G a , so that each 
new value of H a (and thus G a ) corresponds to a new running coupling in the theory. 
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the differential equation (|3.7J1 is inhomogeneous, and, even if g 2 is chosen vanishingly 
small at some scale, the non-zero value of g 2 generates a non- vanishing g 2 below that 
scale. The origin of this is the term proportional to l/£ in Zi ong , cf. Eq. ()3.6j) . 

It turns out that the RG equation (|3.7|) can be solved analytically. This will be 
the topic of the next section. 



4. Running couplings and A parameters 

We begin by rescaling 

t -> 16n 2 t = logO/A) , (4.1) 

in order to get rid of the factors l/(167r 2 ) multiplying both beta functions. The 
solution of Eq. (|3.1|) can then be written as 

where we absorbed the integration constant into the scale A. In other words, we 
define A as the scale where g 2 diverges. We substitute this solution into Eq. ([3.7)1 . 
and introduce a new variable 

x = 9 2 + ^~ 9 2 = 9 2 + 3g 2 ■ (4.3) 
2d 

Equation (|3.7jl then takes the form 

(4 - 4) 

in which 

r 2 be 3^ 

6 = 6 * + * = £tf-3Cfc>0. (4.5) 
4d 2d 2 

The positivity of 6 follows because Ch < AT/2. Equation (|4.4jl has the form of a 
Ricatti equation [TT], and can be exactly solved; its solution is (t > 0) 

x(t) = I Pa+r+ + qa ~ ta ~ , (4.6) 
V ^ dt pt a + + qt a ~ V ; 

where p and q are two arbitrary constants (not both equal to zero), and a± are the 
solutions of a (a — 1) + bd/b 2 = 0, 



1/ / bd\ 1 ±fN 2 -lUNC H + 27C H . 

a * = 2[ 1± r- i v) = 2( 1± V- i urn — £ l- < 4 - 7) 

Obviously, we are looking for real solutions, and this is what we find if p and q are 
chosen real and the radical in Eq. ()4.7|) is non-negative. If the radical is negative, the 
solution of Eq. (|4.4j) may be written differently. Introducing 



W^-4' (4i 



8 



the solution is 

x (t) = -A- ( 1 + - 21 ^ ) , (4.9) 

with an integration constant. 

The two types of solutions, Eqs. (J4.6j) and (J4.9j) are rather different mathemat- 
ically. Taking H to be the Cartan subgroup U{l) N ~ l , for which Ch = 1, one 
finds that for A > 4 the radical is positive, and thus the solution ()4.6|) applies, 
whereas for N = 2 and A = 3, Eq. ()4.9|) is the relevant solution. If we take 
H = SU(M) x SU(N — M) x 17(1), the radical is always negative (with 7 = 
for all iV and M). All solutions are decreasing with increasing t, because dx/dt is 
always negative according to Eq. ()4.4|) . Before we turn to a discussion of each type 
of solution, we note that in both cases, if g 2 = x — 3g 2 is chosen positive at some 
scale t = t , it is guaranteed that g 2 will stay positive at all lower scales down to the 
first value of t where the solution diverges. This is of course because dg 2 /dt is always 
negative. 

4.1. Case 1 - 4 ff > 

When the radical is positive, we have that 1 > a + > \ > a_ > 0, and, with 




x = yi-4-<i, ( 1.10) 

the solution can be written as (t > 0) 

= ^ p(l + ,>v + ,(l~ X ) 
2dt ptx + q 

Depending on the sign of the integration constant q/p, this solution shows two types 
of behavior. When q/p > 0, this solution is always positive for t > 0, and diverges for 
t — *■ 0. In other words, the coupling g 2 diverges at the same scale as the coupling g 2 . 
When q/p < 0, the solution already diverges for a positive value of t with t x = —q/p, 
and the coupling g 2 diverges at a scale A larger than the scale A at which g 2 diverges. 

The sign of q/p can be translated into an inequality between the "initial condi- 
tions" #o = 9 2 (to) an d go = g 2 (to) where to > determines the ultra-violet cutoff 
Ao = A exp(167r 2 to)- In terms of xq = x(to), we have 

_ q = tjjdxpto - a+) ^ 
p dx t — ct_ 

Requiring that — q/p be positive therefore translates into (recall a_ < aq_) 

a + a + b 2 a_ aJb 2 

x > -j- = —j- g , or x < -=- = g . (4.13) 
dt d aiQ d 

As can be seen by combining Eqs. 1)4.11)1 and 1)4.12)1 into 

= j_ dxot (a + t* - q_t^) - a + a_(t x - tg) 1 
dx t (t* - t^) - (a_t^ - a+to) 
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the second of inequalities ()4.13j) leads to a divergence of x(t) at some t > t ; hence the 
only divergence at t < t occurs at t — 0, as for q/p > 0. The first of inequalities (|4.13|) 
leads to a divergence of x(t) at t — t with < t < t , and thus represents the new 
class of solutions with an infra-red scale A = A exp(167r 2 t) such that A < A < A . We 
note that the first of inequalities ()4.13|) implies that certainly xq > 3g$, as required 
by Eq. (jOjl . because a+ > §, and thus a+b/d> 22N/(3N + 6C H ) > 11/3 > 3, using 
< C H < N/2. 

Summarizing, we find solutions diverging at a value of t between zero and to if 
and only if 

^ > -L (2&a+ - c) <? 2 , (4.15) 

and solutions only diverging at t = on the interval < t < to if this inequality is 
not satisfied. One can verify that i {2ba + — c) g 2 > |g 2 . 

4.2. Case 1 - 4 || < 

When the radical is negative, the solutions behave very differently from those dis- 
cussed above. The solutions have an infinite number of zeros, 

x(t' n ) = iff 7log(0 = + £ + n7r , (4.16) 

where (3 = arctan(— 2 / y), and an infinite number of divergences, 

x«) = oo iff 7 log(O = + n7r, (4.17) 

in addition to the divergence at t = 0. Clearly, the zeros and divergences alternate, 
consistent with the fact that x is always decreasing. In terms of the "initial condition" 

Xq = x(t ), (f) can be expressed as 

= 7 log(to) - arctan I 7 . (4.18) 

ydtoxo - 5 y 

The solution x(t) can be written in terms of xq and to by combining Eqs. (|4.9J) 
and JUIHJ): 

x(t) = i ( l - + 7 Jx0t °'r 7tan(7l ° g(t/t0)) ) ■ (4-19) 
dt \2 (dx t - 2) tan (7 log (t/t )) +7/ 

We may now define the infra-red scale A by choosing t as the first value below to 
where x(t) diverges. It follows that always A > A (t is always larger than zero, 
because arbitrarily small positive solutions of Eq. ()4.17|) can be obtained by taking n 
large and negative). 



5. Discussion of implications 

In this section, we discuss the possible implications of the one-loop RG analysis of 
the previous two sections. We emphasize that this discussion is speculative, because a 
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one-loop calculation cannot really tell us what the actual non-perturbative dynamics 
of the theory looks like. 

The first observation is of course that both g 2 and g 2 are asymptotically free. If 
we choose small values for both couplings at some high scale to, they will grow as 
one goes down in energy. The dynamics will be non-perturbative at low scales, while 
perturbation theory should give a reliable description of the dynamics at scales well 
above the larger of A and A. Clearly, in order to find out more, non-perturbative 
techniques, such as numerical lattice computations, will be needed. However, we may 
contemplate what the resummed one-loop calculation, embodied in the solution to 
the RG equations, suggests for the dynamics governed by the coupling g 2 . 

One important question, which we already raised in the Introduction, is whether 
SSB of eBRST symmetry might take place due to the fact that g 2 becomes strong. 
For example, it is not excluded that some ghost condensate is generated, reminiscent 
of the way a chiral condensate is generated in QCD. A question is then at what scale 
this SSB, or any other strong dynamics, would take place. 

The results of our one-loop RG analysis give some interesting hints as to what to 
expect. Let us assume that the message of this analysis is the following: The scale 
at which a running one-loop coupling diverges provides a valid estimate of the scale 
where that coupling becomes strong. The RG analysis thus implies that g 2 cannot 
become strong without g 2 becoming strong too, since the scale A associated with g 2 
is always larger than or (approximately) equal to A, the scale associated with g 2 . 

First, consider the simpler case analyzed in Sec. 4.1. There are two possibilities: 
A = A if #q is chosen small enough relative to g 2 ,; otherwise, A > A (cf. Eq. (|4.15j0 . 
There is no upper limit to A (as long as it is below the cutoff scale A which, in 
turn, can be arbitrarily large). Of course, the precise distinction between these two 
possibilities results from our definition of the infra-red scales as those where the one- 
loop couplings become infinite. Had we defined the infra-red scales by requiring the 
running couplings to take some finite value of order one, we would still find solutions 
with A ~ A, and solutions with A ^> A, but no solutions with A C A. In other words, 
it is likely that there is a region in the phase diagram where the couplings become 
strong simultaneously, and a region where g 2 becomes strong while g 2 is weak, with 
the possibility of choosing A/A arbitrarily large. 

For the case analyzed in Sec. 4.2 the situation looks rather strange at first. If we 
demand that g 2 (t) assumes some fixed finite value for some fixed t > 0, the cutoff 
to cannot be taken arbitrarily large; in fact to must be kept below the first value of 
t where the solution vanishes. In practice, however, the physical interval between 
one value of t where the solution diverges, and the closest higher value where the 
solution vanishes can be very large. The reason is the logarithmic dependence on 
t, which represents double-logarithmic behavior in scale. Two numerical examples 
(taking G/H = SU(2)/U(1)) will suffice to illustrate this. In both examples we take 
A = 100 MeV, and t = 0.29, corresponding to A ~ 10 19 GeV. 

In the first example, we choose xo = l/(2dt ) = 0.86 or gl = 0.16, and it follows 
from Eq. (PLTTIJ) that 

x(t) = X- (1 - 2 7 tan( 7 log(t/t ))) . (5.1) 
2dt 

The first divergence below to happens when 7log(t/to) = — 7r /2, which translates into 
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a value for A about 3% larger than A. So, we can easily arrange that g 2 becomes 
strong at about the same scale as g 2 , while choosing the cutoff at the Planck scale. 

In the second example, we raise x to 1.7 or g 2 = 1.0, which is probably still 
a relatively weak coupling. This leads to A = 75 GeV. We see that, given a cutoff 
scale of 10 19 GeV, and choosing a typical QCD scale for A, we can arrange for A to 
be virtually any scale in between. In practice, the solutions of Sec. 4.2 are not so 
different from the mathematically simpler solutions of Sec. 4.1. 

6. Conclusion 

We calculated the one- loop beta function for the longitudinal coupling g 2 = £g 2 in 
an equivariantly gauge-fixed SU(N) YM theory in which only the Cartan subgroup 
U(1) N ^ or the subgroup H = SU(M) x SU(N - M) x 17(1) is left unfixed. We 
found that also the longitudinal coupling g 2 is asymptotically free, just like the gauge 
coupling g 2 . These results readily generalize to other subgroups, cf. Appendix A. 

We then solved the corresponding one-loop RG equations. Since both couplings 
are asymptotically free, dimensional transmutation takes place for both couplings, 
and an infra-red scale is dynamically generated for each. In the case of the YM 
coupling g 2 we know that the one-loop beta function correctly "predicts" that the 
coupling indeed becomes strong at low energies. Our one-loop results suggest that 
also the longitudinal coupling g 2 becomes strong at low energy when, or even before, 
the gauge coupling g 2 becomes strong. A priori, a possibility might have been that 
g 2 stays weak also when g 2 becomes strong, but this possibility seems excluded by 
the one-loop RG analysis. 

The one-loop RG analysis suggests that there is always some strong-coupling 
dynamics (possibly accompanied by SSB and/or the formation of bound states) in 
the longitudinal sector. As already discussed in the Introduction, a key question is 
whether there are circumstances in which the strong dynamics at the scale A affects 
the gauge-invariant sector. If so, a further question would be how the modified 
dynamics in the gauge-invariant sector depends on the scale A; i.e. whether there 
is any qualitative distinction between the cases A ~ A and A ^> A. Do these two 
possibilities correspond to different phases? This is, to us, a question without an 
obvious answer. In the case that A ^> A, g 2 becomes strong while g 2 stays weak, and 
the strong-coupling dynamics is thus entirely determined by that of the longitudinal 
sector (the "reduced" model, cf. Sec. 1). In the case that A ~ A, all degrees of 
freedom may play an important role, with possibly different implications. 

Clearly, the full dynamics of this gauge-fixed theory, including that of the longi- 
tudinal sector, will have to be further investigated by non-perturbative methods. In 
addition to providing the first hints of what the dynamics of equivariantly gauge-fixed 
YM theories might look like, the one-loop beta-function analysis presented here pro- 
vides a useful first step in that it predicts the scaling behavior in the limit in which 
the cutoff is taken large. 
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Appendix A 

Here we prove Eq. ()3.3j) and compute H a for a large class of subgroups 
H D [/(l)^ -1 , still taking G = SU(N). We will assume that the generators of H 
fit into K square blocks along the main diagonal. The size of the r-th block is 
M r x \ I,., where 1 < M r < N. The number of blocks is constrained by 2 < K < N . 
In the special case K = N all blocks are lxl, which corresponds to H = [/(l) 1 ^^ 1 
P]. In this appendix we thus assume K < N, implying that there is at least one block 
with size M r > 1, corresponding to a subgroup SU(M r ). The special case K = 2 
corresponds to the maximal subgroup SU (M) x SU(N — M) x £7(1) discussed in this 
paper. The case K — 1 is excluded as it corresponds to H = SU(N) = G. Note that 
the number of Cartan generators of SU(N) contained in the non-abelian subgroups 
(the blocks with M r > 1) is Y^=\{M r — 1). It will be convenient to assume that the 
blocks are ordered by increasing size: M\ < M2 ■ ■ ■ < Mk, implying in particular 
that any lxl blocks occur before all the M r > 1 blocks. 

We next construct a convenient basis for the off-diagonal SU(N) generators, 
which includes all coset generators. We define N x N matrices r\ B where 1 < A < 
B < N and 1 < k < 3. The matrix r\ B has precisely two non-zero entries, each with 
absolute value equal to 1/2, chosen such that t\ b reduces to (half) the corresponding 
Pauli matrix, <J k /2, if we keep only the A-th and 5-th rows and columns. The coset 
generators are spanned by t% b , where s = 1,2, and where the non-zero entries are 
constrained to be outside of all the diagonal blocks. For the maximal subgroups 
(assuming M < N/2 to conform to the above ordering), this implies 1 < A < M and 
M + 1 < B < N, and it is easy to check by inspection that f a p^ = 0. In this case, 
one thus finds G a = Cq = 0. 

It is convenient to make a unitary transformation to the non-hermitian rj[ B = 
( T AB { t ab)\ which satisfy the orthogonality relation tr (rJ B Tq D ) — 

\5ac$bd- The matrices r\ B each have a single non-zero entry located above the main 
diagonal. Using this basis, Eq. (J3.3)) may be re-written as 

- 2 X) tr ([T* , rt B ] [T\ r£,]) = H AB 8 AC 5 BD . (A.l) 

i 

Denote the number of 1 x 1 blocks by L. (Our assumptions imply < L < N — 2.) 
Eq. (jA.l|) is proved by considering the following three cases: 

• 1 < A < B < L (this requires L > 2). This effectively reduces to the case of 
H = [/(l)^ -1 . Constancy of Hab follows from skewed permutation symmetry 
PP. The actual numerical value is found to be H A b — !■ 

• 1 < A < L < B. We say that B "belongs" to the r-th block if the 5-th entry 
along the main diagonal resides in the r-th block. By assumption, the size of this 
block is M r > 1. It is easy to see that t\ b transforms as the anti-fundamental 
representation under SU(M r ). We calculate Hab by summing the contribution 
from all the Cartan generators of SU(N), which is always equal to one, plus the 
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contribution of the extra, off-diagonal SU(M r ) generators. We use the known 
value of the quadratic Casimir of the fundamental representation, as well as 
J2iTi = (N — 1)/(2N), which is true for the sum over the Cartan generators 
of any SU(N) in the fundamental representation. This proves Eq. (|A.1|) . with 
H AB = 1 + D(M r ) where D{M) = (M — l)/2. 

• L < A < B. Now A and B belong to the r-th and s-th blocks respectively, where 
1 < M r < M s . The matrix r AB transforms as the fundamental representation 
of SU(M r ), and as the anti-fundamental representation under SU(M S ), and 
Eq. ipQjl follows with H AB = 1 + D(M r ) + D(M S ). 

As a check, consider H = SU(M) x SU(N-M) x U(l). For M > 2, this corresponds 
to the third case, and H AB = 1 + D(M) + D(N - M) = N/2. For M — 1 one has 
H = SU(N — 1) x £7(1), which corresponds to the second case, and again H AB = 
1 + D(N — 1) = N/2. Another example of a subgroup with a constant H AB , which 
is a generalization of the case H = ^(l) 7 ^ 1 , is H = SU{M) K x [/(l)^ 1 with 
G = SU(N = MK). In this case, the subgroup's generators fit into K blocks of size 
M x M, and H AB = M. 
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